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Abstract 

Using modern data, we first refit constants needed to complete an anomalous irpuai 
Lagrangian obtained within the approach of hidden local symmetries. Then we de- 
rive from this Lagrangian electromagnetic isoscalar pir^f and paij exchange currents 
needed in calculations of the deuteron electromagnetic form factors at large momen- 
tum transfers. 

Key words: electromagnetic, isoscalar, exchange, current, anomalous, chiral, 
Lagrangian 

PACS: 11.40.-q, 12.39.Fe, 25.30.Bf 



1 Introduction 



In calculations of the electromagnetic form factors of the deuteron, the prey 
exchange current plays an important role [1-3]. According to [4-6], it is derived 
from a prey vertex 

<n m (q 2 )\jr-\p l M)>= j 6 WWg 2 ) ^ Mm; , (1) 

" m p 

where q = q 2 — q\ and the constant g pvr7 is extracted from the data on the 
width of the decay p — > n + 7. 
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Here we consider the task of constructing this current by using an anoma- 
lous Lagrangian. Such a Lagrangian of the irpuaifi system was constructed 
in [7] within the approach of hidden local symmetries . It includes also the 
external electromagnetic field. The subsystem npuai was considered later in 
[8], with both the external electromagnetic and weak vector and axial-vector 
fields included, however. In the Lagrangian, several terms are present, which 
yield, besides the current (1), a correction to it and also a new current pai'j. 
The Lagrangian is characterized by several parameters, which should be de- 
termined from the data on various reactions. Experimental data available 10 
years ago allowed Kaiser and Meissner [7] to extract only some of them. The 
present experimental situation [9] and results of the recent work [10] allowed us 
to improve the analysis. As a result, all 4 parameters entering the Lagrangian 
of the system npuai are now available. 

We use further the constructed anomalous Lagrangian to derive the electro- 
magnetic isoscalar pixj and pa\j exchange currents, which can be employed 
in calculations of the electromagnetic form factors of the deuteron. They dif- 
fer from the standard approach by the presence of an additional momentum 
dependence in the electromagnetic form factor of the pivy current and by the 
completely new pa\j current. We would like to stress that the considered 
model fully respects chiral invariance and vector dominance and it is consis- 
tent with the present experimental knowledge of elementary processes such as 
radiative decays p — > nj, uj — > ttj, f\ — > pj and f\ — > pixix. 

In Sect. 2, we introduce the anomalous Lagrangian, extract the needed con- 
stants from the data and present the electromagnetic form factors in our model 
for the vertices pirj and pa^. In Sect. 3, we construct the associated electro- 
magnetic isoscalar exchange currents and make the non-relativistic reduction 
of these currents. The discussion of the results and our conclusions are given 
in Sect. 4. 



2 Anomalous Lagrangian of the npuoai system 

We derive the pnj and pa\j vertices and the associated exchange current from 
an anomalous irpujcii Lagrangian [7,8]. Besides the meson fields, the external 
vector isoscalar £> M and isovector V M and axial vector isovector fields are 
included. The Lagrangian, which is of interest for the problem under investi- 
gation, reads 

C 7 = 2ige KX ^ {d K u x [{gpp - eV M ) • d u n + eA u )] 
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+ (gu K - \eB K ) [(d x p n ) • 9,7? + eX)]} , (2) 
£ 8 = -2ige K \^ {d K uj\ \(gp p - eV M ) ■ (#4 + -^-<9„7r)] 

1 1 

+ (^k - ^eB K ) [(dxPp) • {ga v + — d v W)\} , (3) 
C 9 = 2iee K x^ u {\d K B x \{gp^ - eV M ) ■ (— d v ir + eX)] 

+ (^ K - JeB K ) [(9 A V M )(^^7f + eX)]}, (4) 
£io = -2ie£ K A/^ { ^<9«#A - eV /t ) ■ (ga v + ^^tt)] 

+ (gu K - ^e£ K ) [(9 A V M ) • (gd v + ^d u n)]} . (5) 

Here g = g p and e is the elementary charge. It follows from the KSFR relation, 
2/2 g 2 = m 2 ; that the yalue f g is ~ 6, if / w = 92 MeV and m p = 770 MeV 

are used. 

The strong vertices poo 7r and puj a\ are already known from the appendix A of 
Ref. [7]. Then Eqs.(2)-(5) provide all other terms arising from the homogenous 
terms of the anomalous action due to the presence of the external electroweak 
interactions which change the natural parity. 

The total Yukawa-type anomalous Lagrangian of the npuja\ system is given 
by the sum 

Can ^ ] Q C-i ; (6) 

i=7,10 



where the constants, q, according to [8] are 

C 7 — C 7 + -Cg , Cg — Cg , Cg — Cg + —Ciq , C\q — C\q , (J) 



and the constants q are given in [7] as Q 
g 7 = 1.32 x 1(T 2 , c 8 = -2.05 x 10" 1 , c 9 = -5.14 x 10" 3 , c w = 



However, new data on several reactions relevant to the analysis aiming to get 
these constants has recently been published [9] and also a new work [10], in 
addition to [7], has appeared, which enables us to get more reliable values of 



1 The constants 5j used here differ from those of Ref. [7] by the factor g. 
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them. But we still cannot remove an uncertainty in the sign of the constants, 
implicitly present also in the previous studies [7,10]. 



2. 1 Extracting the constants from the data 



We start by extracting the constant c 7 . It enters the effective constant g pu}7T 
defined as 

£pwir = i^f^^KXpudnUxPp- d u if . (9) 

Jtt 



The constant g puJ7T was found [10] to have the value 

^ = 1.2. (10) 

On the other hand, our model Lagrangian (6) leads to 

9 pun = 4g 2 c 7 , (11) 

which yields 

c 7 = 8.64 x 10" 3 . (12) 

This value of c 7 is lower by ~ 30% than its older value given in (8). 

Using the Lagrangian (6) we derive the amplitudes for the radiative processes 
p — > 7T7 and uj — > with the constant Cg entering in combination with c 7 
into the effective decay constants 

2gm p _ 2#m 
dp*-/ = ^r^{c 7 + c 9 ) , g unl = — — (c 7 + c 9 ). (13) 

These decay constants enter the radiative decay width of a vector meson V 
decay into a pseudoscalar meson P as [11] 



T-i a 2 

1 v ^ p 7 - ^g VPl m v 



m P \ 2 



1- — 

\mvJ 



(14) 



According to [9], the widths are 

rp±_ w ± 7 = (67.82 ± 7.54) fceF , r w _^ 7r o 7 = (714.9 ± 42.1) keV . (15) 
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From Eqs. (12)— (15) we have 

c% = 8.64 x 1(T 3 , Cg 1 = 9.55 x 1(T 3 . (16) 

The superscript c (n) implies that the corresponding quantity is obtained from 
the charged (neutral) meson decay. The values Cg and c"g differ by ~ 10%. 

We calculate the weighted value of Cg according to the equation 



a w = \ Y, u>i / E w * I > w * = r */ Ar < • ( 17 ) 

\i=l,n I \i=l,n 



Using Eqs. (15) and (16) we obtain 

= 9.23 x 10" 3 . (18) 

In comparison with eg from Eq. (8), the new value (18) has the opposite sign 
and it is larger by a factor of ~ 2. 

If we admit negative values for the sum £7 + c 9 then we get 

Cg° = -2.65 x 1(T 2 . (19) 

With the weighted value (18) we get for the radiative decay widths 

rp±_ w ± 7 = 72.5 keV, r w _ ff o 7 = 689.6 keV , (20) 

and for their ratio 

r a ;^ 7 ro 7 /r p ±^ 7r ± 7 = 9.5, (21) 

and for the effective decay constants 

g pwi = 0.585, ^ = 1.78. (22) 

For the negative values, the effective couplings have the opposite sign. 

We can also use the same procedure with the neutral p meson decay, for which 
the decay width is [9] 

r p o_^o 7 = (102.5 ± 25.6) keV. (23) 



The result for the constant eg from this equation is 

a 9 = 1.25 x 1CT 2 . 



(24) 



The experimental value of r p o_ +7r o 7 is converging in time to the value of 
r p ±^ 7r ± 7 [9], but it is still impossible to get both decay widths at the la 
error level by using the same value of eg, since the error bounds do not over- 
lap. We can use Eq. (17) for the neutral meson decays to get the weighted 
value of Cg = 1.014 x 10~ 2 , but the calculated value of r w ^ +7r o 7 = 762 keV 
is already outside the la error bound 757 keV. It can be verified that only 
the values of c 9 from the interval 0.97 x 10~ 2 < c 9 < 1.05 x 10 -2 yield 
simultaneously acceptable values of r u _> n o 7 and r p o^ 7r o 7 . For 

g 9 = 1.00 x 10~ 2 , (25) 

one gets 

= 751 keV, Tpo^^ = 79 A keV , (26) 

with the ratio 

r^^^Oj/FpO^^ = 9.5, (27) 

and for the effective decay constants a value of 

g p7ry = 0.610, <Wy = l-86. (28) 

One obtains similar results with the negative value of c 9 

c 9 = -2.73 x 10~ 2 . (29) 



Actually, the decay widths r^^ 7r o 7 and r p o^ 7r o 7 were used in [10] to extract 
the constant d, related to our constants by the equation d = g (c 7 + c 9 ), with 
the result d~0.01. On the other hand, taking our constants g, c-j (12) and 
eg (25) we get dw0.ll. Evidently, the calculations [10] give the value of d 2 
instead of d. 

It is interesting to follow the change of the value of the constant g W7r7 in 
time. The first use of the pn'j exchange current in the calculations of the 
electromagnetic deuteron form factors and an estimate of the constant g W7r7 
was made by Chemtob, Moniz and Rho (CMR) in [5]. Actually, CMR adopted 
the pnj exchange current constructed by Adler [4]. The existing data on 
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[12] allowed CMR to extract the upper bound | g p7r7 | < 1.02. Inspite of an 
estimate, based on the existing relativistic quark model [13] , that the value of 
g p7r7 migt be about 50% of the upper bound, CMR adopted the value g pvr7 = 1. 
Afterwards, the value g p7r7 = 0.34-0.48 was extracted by Gari and Hyuga [6] 
from the more precise data [14]. Later on, Towner [15] analysed the data [16] 
to obtain the value g p7r7 = 0.578, which is larger by about 50%. At last, in 
[18,19], the value of g p7r7 , obtained from the data [20], has oscillated down to 
0.56. However, our value (22) of g p7r7 = 0.585 is again larger and close to the 
one obtained in [15]. Towner also extracted from the data [17] the constant 
g^7r 7 = 1-98, which is larger by about 10% than our = 1.78 (22). 

Let us note that our analysis is based on the data recommended recently in 
[9] and that in our approach, we extract g p7r7 and g W7r7 from the data simul- 
taneously. This is due to the fact that in our model, both vertices depend on 
the same parameters c 7 and c 9 . 

We conclude that the present experimental results allow us to extract two 
values of the constant eg which differ by ~ 10% 0. However, the difference in 
the effective coupling constants g p7r7 and g^^ given in Eqs. (22) and (28) is 
only 4%. One can also ignore the data on the radiative decay p° — > 7r° + 7 
and use only eg from Eq. (18) and the decay constants (22) with the hope that 
the experimental study of this decay will improve. 

Let us turn to determining the constants c 8 and c w , needed to calculate the 
anomalous processes with the ai meson. However, the decays of this meson 
are not well measured and we use the fi meson decays to extract them from 
the data. In contrast to the experimental situation in late eighties, the decay 
/1 — > pmr is now well measured [9] and good data for the radiative decay 
fi — > p7 has appeared [9]. According to [7], one of the contact couplings, 
h'f ipnn , is connected to c 8 as 

h' 

C8 — 7) — 

g 2 

Using Eq. (4.20) of [7] with the new value 

T fl ^ pwn = (2.808 ± 0.360) MeV , (31) 

one gets 

~ti fif)7T7T = 3.775, (32) 



2 In the case of the negative values, the uncertainty in eg is about 3%. 
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(30) 



and from (30) 



c 8 = 



-1.02 x 10" 1 . 



(33) 



Again, the new value of c 8 is lower by a factor of ~ 2. 

In order to calculate the radiative f\ — > p7 decay, we adopt the pertinent 
part of the anomalous Lagrangian (3.12) of Ref. [7], 



A£ h = cgA£ 9 + ci A£i2, 

A£ 9 = -2ig 2 e KXf3l/ f lK [g(dxPp) • p v - {d x p Q ){eB v )\ , 
£12 = -2ig 2 e K x(3ufi K (d x eB^)p° u . 



(34) 
(35) 
(36) 



With these vertices, we get for the radiative decay width an equation similar 
to Eq. (14) 



a 



1 + 



m, 



(37) 



where 



9phi = 2 9 ( c 8 + Cio) 



(38) 



Using the experimental value of r^-^ [9] 

Tf^n = (1.296 ± 0.288) MeV , 



(39) 



we get 



g 8 + cio = 2.65 x 10~ 2 , 



(40) 



and for the constant cio the value 

cio = 1-29 x 10" 1 , 



(41) 



which is in absolute value even larger than c 8 , Eq. (33). Then for the effective 
coupling <7/ iP7 we have the value 



2/1P7 = L84 - 



(42) 
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Admitting negative values for the sum eg + cio we get 

c 10 = 7.59 x 1(T 2 . 



(43) 



We now have in Eqs. (12), (18), (25), (33) and (41) all the needed constants 
and we can now proceed to constructing the prrj and pai7 exchange currents. 
However, let us first discuss the electromagnetic form factors provided by our 
model. 

2.2 Electromagnetic form factors 

The Lagrangian of Eq. (6) yields the vertex pn^, Eq. (1), with the form factor 

\ Cg -\- Cj / 



1 - F^q 2 ) 



(44) 



where the form factor F^q 2 ) is given by the vector dominance. Since eg > £7, 
the second term on the right hand side of Eq. (44) gives a 3-7% positive 
correction to the first term at q 2 ~ m 2 . In the case of the negative values of 
eg, this term is larger by a factor of ~ 2 than the first one at these values of 



-»2 

Q ■ 



The presence of the terms Cs and £10 in Eq. (6) shows the existence of the 
vertices pai7, which give rise to the corresponding pai7 currents. In analogy 
with Eq. (1) we have 



< <u(l2)\Jx m -\pl(qi) >= eg paiJ e Xl/af , [K^q^q. + K 2 (q 2 )q la \ 5 ml , (45) 



where the form factors Ki :2 (q 2 ) are defined as 



K^q 2 ) = -: 



C8 + Cio 



cs-FU"? 2 ) + cio , K 2 {q 



c 8 



c 8 + Cio 



F w {q 2 ) - 1 . (46) 



The effective coupling constant g pail is 



(47) 



with the constant <7 P / 17 given in Eq. (38). This gives a new contribution to the 
current. 

In the next section, we derive the pn^ and pai7 exchange currents. 
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3 The p7T7 and pa^ exchange currents 



The general structure of the exchange currents which follows from our La- 
grangian (6) is given in Fig. 1. The currents of the pion range are, 



JT, C = ^9,NN^ 1 e, K(3l/ F UJ (q 2 )q K V^l,2)q 2l/ + (1 <- 2) , (48) 
J Tt fi = -i9«NN g2{C l~ ° g) £^uq K Vp{l,2)q 2v + (1 <- 2) , (49) 



where 



7^(1,2) = r? i/3 A^( ?1 2 )A-( 9 |)r« i5 , (50) 



and 



K T 7 

r ?,/j = «(p|)(7/3 - 7^V(36Ks)T a u(pi) and r« 5 = utf^T'uipj) . (51) 



The sum of the currents given by Eqs. (48) and (49) is 

J*p = ig-KNNT— 1 K p7ry (q 2 ) e pK/3u q K Vp{l,2)q 2u + (1 <-> 2) , (52) 

Z Tflp 



where the effective coupling constant g pnl and the form factor K p7T ^(q 2 ) are 
defined in Eqs. (13) and (44), respectively. 



Pi 



Pi 



Bi 

qi 



a 



P 2 



P2 



Pi 



Pi 




P2 



P2 



Fig. 1. The general structure of the electromagnetic isoscalar pTrj and pai7 
exchange currents J^q) considered in this paper. The meson B\ is the p meson. 
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The range of the current is given by the meson B 2 which is either it or a\ 
meson. The graph a is for a mesonic current J^ai where the electromagnetic 
isoscalar interaction is mediated by the meson B which is uo meson. Contact 
terms are given by the graph b, J% 2p1 where the isoscalar electromagnetic 
current interacts directly with the mesons Bi and B 2 . 

The currents of the a± meson range are 



JZ, C - = ~9%e^ v F u ^)(q K + q lK )V^(l,2) + (1 <- 2) , (53) 



= 9 -fg A e^u (csQi. ~ ~doq K ) Vp v (l,2) + (1 <- 2) , (54) 
where the weak axial coupling constant qa = 1.26, and 

7%(1,2) = Yl^^ F {ql)^{qi)T a 2 ^andTl bl/ = u{p' 3 ) lul ,r a u{p 3 ) . (55) 

The total pai'j current is 
Ja^ = -^g 2 g P a 11 £»Kp v [i^i(g 2 )g K + K 2 (q 2 )q lK ] Vp v {l,2) + (1 <-> 2) , (56) 

with the effective coupling constant g pail and with the form factors K 1>2 (q 2 ) 
given in Eqs. (47) and (46), respectively. 

In deriving Eqs. (48) and (49) and Eqs. (53) and (54), we use the equation 
g u = 3g p , which ensures the current conservation. 



3.1 Non-relativistic reduction 



Now we proceed to the non-relativistic reduction of the currents given by 
Eqs. (52) and (56). Having in mind the application of our results to the elastic 
electron scattering by nuclei, we neglect the energy transfer putting q = 0. 
This reduction yields for the space component of the currents 



X = CV ( (<f x q 2 ) 



4M 2 



^ + kv ) (q x q 2 ) <?i 2 + i °\ • Pi x qi 



-1 (1 + K V ) (q 2 ■ P 2 ) (qx {3 X X £))}) 
• (* 2 ■ q 2 ) (n ■ r 2 ) A F (q 2 ) A F (q 2 ) + (1 



2), 



4 = i g -^g 2 g pail (K 1 (q 2 )q + K 2 {q 2 )qi) x (b 2 {\ 1 2/M 



8M 2 



-> 2 
0i 



+i a 1 -P 1 xq 1 ]}-i (ai x y (<r 2 • P 2 ) 



AM 2 



(57) 
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(fW^A^A^I) + (1 - 2). 



(58) 



Here 



Ctt — — % 



Qttnn g 



2M 2m, 



(59) 



and 



Pi = Pi + Pi- 



rn 



The importance of relativistic corrections in the space component of the pii'-f 
current was stressed in [18]. They were obtained earlier in [21]. We present 
them also in the a\ exchange term. An evaluation of the contribution of these 
terms to the electromagnetic deuteron form factors would give more confidence 
in the validity of the static approximation for the current (58). 

Besides the relativistic corrections present in the current (57), we list the boost 
current J® arising from the commutator of the leading terms in Eq. (57) with 
the operator of the kinematical boost and the retardation current J^ et . The 
kinematical boost current is defined as [22] 



J, 



B 



(p) , Jl 



(61) 



where 



Xo(P) = 



16 M 2 

f = r i - r 2 , 



-P 2 - 2i (f • P) (p- P) + 2i (<ri - <r 2 ) x p ■ P] , (62) 

(63) 



and the static approximation to the pn^f current J® is defined as 

JJ = a (<? X q 2 ) (B 2 ■ q 2 ) (n ■ r 2 ) A p F (q?) ^M) + (1 - 2) . (64) 

Generally, the current is quite complicated. However, it simplifies in the 
Breit system, where it reads in the cm. coordinates, 



Ji 



■C n (q x q_) \ (a 2 ■ q 



16M 2 

- • K x qj + i (a 2 x q_) ■ (k x qj } 
• (n ■ r 2 ) Ap (ql) A F (f_) + (1 - 2) . 



-j(r-q) (q-k) + (ffi ■ K x qj 



(65) 



Here 
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\ (Pi - #>) , P = Pi + P2 , P' = Pi ~ P2 , 



k = j? -p, K = j? + p, q = P'-P, Q = P' + P. 
Qi = Pi - Pi = k + -q = q+ , & = fi 2 -P2 = -k + ^q 

Pi=ti+Pi = K+^Q, P2 = 02+P2 = ~K+ l -Q. 



-Q- 



(66) 



In the Breit system, Q = P' - P = and P' 



q/2. 



In the realistic calculations, the currents, Eq. (52) and Eq. (56), are multiplied 
by the strong form factors. In essence, this procedure reduces in the case of 
the p7T7 current to the following change in Eq. (50), 

K 0) A£ {qi ) - F pNN 0) A P F (0) F nNN {0) A F [0) , (67) 



and to the change 

K (0) ^ [0) 



F pNN {0) A P F 0) F aiNN 0) Ap 1 0) , (68) 



in Eq. (55) in the case of the pai7 current. While in this case, the static 
approximation is expected to be reasonably good, for the puj current the 
retardation current arises by the following change in the current J®, 



A£ 0) A F 0) 



F pNN 0) A F 0) F wNN 0) A F yqi 

\K0)-^ln(F pNN 0)) 



x 1 + qf Q 



+Q20 



d 



A F (0)-—ln (F nNN 



(69) 



The retardation current J^ et is proportional to the energy transfers q±o and 
q2o defined as 



2M 



(70) 



which can be expressed via the cm. variables by using Eqs. (66). The part 
of the current J£ et without the contribution from the derivative of the form 
factors was derived earlier in [21]. 



The exchange charge densities are 
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J-rrO — 



J a! — 



Q-kNN 9 1 + Ky 

2M 2m p 2M 



9a 2 i_ 

2 9 9pai1 2M 



Ky 



9 P ^K piT1 (q 2 ) [q- q 2 x (a l x q t )] 
■ (fi • f 2 ) A p F (q 2 )A F (q 2 ) + (1 <-> 2), 
Xi(<7 2 )<7 + K 2 (fqi] • [<? 2 X (ft x ^)] 

• (f^A^A^ 2 ) + (1 - 2). 



(71) 



(72) 



4 Results and conclusions 



In this paper, we have studied the electromagnetic processes in the npaioo 
system, characterized by the change of the natural parity of the participating 
particles, by using the anomalous chiral Lagrangian [7,8]. This Lagrangian 
contains 4 constants c 7 -c 10 , which were first extracted from the data in [7]. 
Using new data [9] on the radiative decays of the p-, u>- and fi mesons and on 
the decay f\ — > pirn and the constant c 7 from [10], we have determined the 
constants cs, eg and cio- The new constants differ considerably from the old 
ones. These constants enter the decay rates in the combinations c-j + cg and cs + 
Cio, which can be obtained only up to the sign. We considered both possibilities 
and we extracted two sets of constants. Disregarding the sign problem, these 
constants are now known within a 10-15% accuracy. The main source of the 
uncertainty is due to an inconsistency in the data for the radiative decays of 
the charged and neutral p mesons and of the u meson. However, the effective 
couplings g fml and g W7r7 suffer from this inconsistency by an uncertainty only 
of about 4%. 

We have subsequently constructed the electromagnetic form factor K pnl (q 2 ) 
that contains, besides the well known term, an additional piece which can be 
large, if the corresponding choice of the constant eg is made. The presence of 
the pa\UJ vertex in the Lagrangian gives rise to the new pai^ vertex, which 
contains two form factors Ki(q 2 ) and K 2 (q 2 ). 

The isoscalar exchange currents, associated with these vertices, are constructed 
as well. In addition to the known pnj current, we have also pa\^ current with 
well defined parameters extracted from the data. This is in contrast to ooej 
[18] or [19] currents based on less firm ground. After the non-relativistic 
reduction, relativistic corrections are also preserved in the space component 
of our currents, in conformity with the results of Refs. [18], [21]. As to the pii'-f 
current, the leading terms on the right hand sides of Eqs. (71) and (57) are 
the same as the exchange charge density (3. 2a) and the space component 
(3.2b) in [6], respectively. However, our electromagnetic form factor K pnl (q 2 ) 
Eq. (44) exhibits an additional momentum dependence. Our relativistic correc- 
tions contain the contributions from the vertices, from the kinematical boost 
and from the retardation, with the part arising from the strong form factors 
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included. The current J ai Eq. (58) and the density J ai0 Eq. (72)are entirely 
new. The subsequent numerical work should reveal how much they can influ- 
ence the deuteron form factors at high momentum transfers. 
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